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INTRODUCTION 

a^ ■ 

Q^ ' Nearly thirty years of chaos research have shown that chaos is an ubiquitous phenomenon. 

^ ■ Almost any dynamical system with a sufficiently large phase space described by nonlinear equa- 

■ tions shows some kind of complexity in its time evolution, i.e. chaos. The degree of complexity 
^ ■ is measured by topological entropy, Lyapunov exponents, decay of correlations, and, in the case 
CN ■ of dissipative systems, various fractal dimensions. Examples abound, ranging from nonlinear 

oscillators to hydrodynamical systems and the motion of planets and asteroids (Schuster 1988). 
Nevertheless, these investigations have also revealed significant differences between an irregular 
looking chaotic signal and a really and truly random one, due to 'order within chaos' (the title 
of a conference, Campbell and Rose 1983, and a book, Berge, Pomeau and Vidal 1988). One 
cr^ ■ particular example of this order and coherence in chaos will be exploited below. 

o : 

m , 

G^ ■ The fascination excerted by 'quantum chaos' and the controversy that surrounds it derives 

■ partly from the fact that according to the above measures, there is no chaos! For once, there 
y^*- is no nonlinearity in the Schrodinger equation. Secondly, quantal spectra of bounded systems 

Q ■ are discrete, so that the time evolution of any observable is multiply periodic. On the other 
^ . hand, the classical time evolution of chaotic systems is described by a non degenerate continuous 
O ■ spectrum, which is capable of rapid dephasing and decay (Cornfeld, Fomin and Sinai 1982). Since 
J> ■ according to Ehrenfest's theorem, classical and quantum observables evolve alike for short times, 
^ . one expects to see differences between classical and quantum behaviour after times T sufficiently 
long to resolve the discreteness of the quantum spectrum. [] One now could try to develop 
measures that characterize the transient chaos-like properties of quantum systems whose classical 
counterpart shows chaos or focus on generic properties of quantum systems in analogy to chaos 
being the generic property of classical systems (e.g. Eckhardt 1988a). 
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Alternatively, one may use the association 'chaos = irregularity' and ask whether one can 
identify within wave theory properties that give rise to irregularities of some kind. A hint is 
provided by wave optics: a superposition of a sufficient number of rays or waves leads to irregular 
speckle patterns. This analogy has been studied by O'Connor, Gehlen and Heller (1987), who 
also found a difference between patterns due to interference of waves with random amplitudes and 

^This is an optimist's estimate. A pessimist might estimate the tm'nover time from the evolution of Gaussian 
wavepackets and conclude that it goes like log(l/?i) because of exponential spreading (Berry and Balazs 1979). 
Recent work by Tomsovic and Heller (1991) seems to support the optimist's version. 

^Since there is no intrinsic definition, the correspondence to a classical system is required, with all its limitations 
and pecularities as for example in the case of spin systems (Graham and Hohnerbach 1984). 
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phases (as in the usual optical speckles) or due to waves with random phases but fixed amplitudes 
(the case relevant for quantum waves in billiards). This line of thought leads to random matrix 
theory (Haake 1991). 

For quantum systems, the equivalent of light rays is provided by the semiclassical propagator, 
which expresses the transition amplitude from one point to another in a given time interval as a 
sum over all classically allowed paths connecting the two points in that time (Miller 1974). This 
semiclassical propagator in the time domain has attracted some attention recently (Tomsovic and 
Heller 1991). In most cases, however, one transforms from the time domain to the energy domain 
and takes the trace which then allows one to study the spectrum of energy eigenvalues of the 
system (Gutzwiller 1967, 1969, 1970, 1971, 1990). In one degree of freedom one ends up with the 
WKB method (for phase shifts, propagators, quantization rules etc.), which is known to provide 
accurate and useful estimates (Berry and Mount 1972). 

In many degree of freedom chaotic systems the semiclassical analysis leads to a relation between 
the quantum eigenvalues and classical periodic orbits (Gutzwiller 1990). The analysis is not easy 
to carry through to the end, with a quantization rule as simple as WKB as the final answer, 
since the number of periodic orbits grows exponentially with period, thus prohibiting a direct 
evaluation (Eckhardt and Aurell 1989). To some extend this exponential proliferation is the 
cause for irregular features in quantum systems, especially in the wavefunctions (Voros 1976, 
Berry 1977). 

However, one also has a semiclassical expression, which in principle should be quantitatively 
useful when Planck's constant is sufficiently small, the de Broglie wavelength sufficiently short. 
The problem is how to tame the exponential proliferation of terms, that is, how to overcome the 
'topological pressure' (Ruelle 1978, Gaspard 1992) or the 'entropy barrier' (Berry and Keating 
1990). This problem occurs in a variety of contexts also outside semiclassical mechanics and 
techniques have been developed to overcome exactly this proliferation problem (Artuso, Aurell 
and Cvitanovic 1990a,b). Below I will demonstrate both the taming of exponential proliferation 
of classical orbits and the emergence of irregular features for the time delay in scattering. I 
will do so for a specific model, scattering off three disks, that has proven extremely useful in 
developping and testing a number of concepts in chaotic scattering (Eckhardt 1987, Gaspard and 
Rice 1989a,b,c, Gvitanovic and Eckhardt 1989, Eckhardt et al 1992, Eckhardt and Russberg 1992, 
Tanner et al 1991). 

WIGNER'S TIME DELAY 

The quantity for which I shall discuss the semiclassical expansion is Wigners time delay (Wigner 
1955). As discussed by Wigner (single scattering channel) and Smith (1960) (many channels), it 
describes the phase difference between a scattered wave and a freely propagating one (though 
with the same asymptotic motion and thus with perhaps a single reflection in the interaction 
region, cf. Narnhofcr and Thirring 1981). It controUs the absorption in a scattering experiment 
(Doron and Smilansky 1992, Doron, Smilansky and Frenkel 1991) and is of some relevance to 
persistent currents in mesoscopic rings (Akkermans et al 1991). 
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In terms of the ^'-matrix, the time delay t(/c) is given by 



dS d 

r ^ -ihtr S^— = -ih—\ndetS . (1) 
oE oE 

The semiclassical expression is identical to Gutzwillers (1990) density of states (i.e. the excess 
density of states above a uniform background) and has been mentioned by Balian and Bloch 
(1974). For the billiard problem studied here, energy and Planck's constant only appear in the 
combination of the wavenumber, k — \/2mE/h (m is the mass of the particle). When expressed 
in terms of wavenumber rather than energy, the time-delay becomes a 'length'-delay, viz. 

oo piLpki — i/ip7rr/2 

p r=i ^|det(l-M;)| 

As usual in periodic orbit theory (Gutzwiller 1990, Eckhardt 1992), the semiclassical expression 
splits into a smooth part (to(/c)) and a fluctuating part (t/(A;)) determined by periodic orbits. 
Since we have expressed everything in terms of wavenumber k rather than energy E, the time 
delay contains the geometrical length Lp of the trapped trajectory, the Maslov phase /ip and the 
linearization perpendicular to the orbit Mp. The summation on r accounts for multiple traversals 
of an orbit. This expression is very similar to the classical trapping time distribution, where the 
amplitude of a trapped periodic orbit would be 1/| det(l — Mp)\, i.e. without the square root and 
phases (Kadanoff and Tang 1984). This difference reflects the one between a quantum amplitude 
and a classical probabiUty. 

For the two degree of freedom systems considered below, Mp has eigenvalues Ap and 1/Ap (by 
convention, |Ap| > 1). Using an expansion of Miller (1975) for the determinant, one can express 
of the time delay as a logarithmic derivative of an infinite product: 

Tf{k) = Re EE— ^ / (3) 

Pi oo oo ^ 1 s 

= Re E E E (-d {e^'^'-'''^^''\K\-"'K'Y (4) 

<^l^ p j=Or=l ^ ^ 

= Re^^EElog(l-e^^'''-^'^-^/'|AprVV') (5) 

OK p 

where the prime denotes a derivative with respect to wavenumber k and 

oo 

= n n (l - e^^^^^-^'^-^/'IApl-i/^^;^-! 
p j=o 
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Fig. 1: Proliferation of periodic orbits with period for the three disk system for 
different distance to radius separations. From top to bottom: d : R = 2.5, 2.8, 3, 4, 5 
and 6. 

is a Selberg type zeta function (Selberg 1956, Voros 1988), named after similar expressions in the 
theory of motion on surfaces of constant negative curvature (Balazs and Voros 1986). 

The technical difficulty with this expression is the exponential proliferation of periodic orbits, as 
mentioned before. For the three disk billiard, this is easily seen from the topology of trajectories: 
After leaving a disk, the particle has two choices for the next collision. Thus the number of orbits 
grows like 3 ■ 2"^^ with n the number of collisions (there are three choices for the first and two 
for the following ones) (Eckhardt 1987, 1992). Taking into account an average length of paths 
between collisions, this translates into an exponential proliferation with length (Fig. 1). 

For practical calculations it is now useful to group terms according to their topological length. 
Technically, one introduces a counting variable z"'p where Up is the topological period of the orbit 
and expands 



in a power series using formulas of Plemelj (1909) and Smithies (1941) (see also Reed and Simon 
1983). In Figure 2 the coefficients Cn{k) for various distance to radius ratios within the Ai 
symmetry are shown. One clearly notes the rapid decay of coefficients with increasing n. 

As a careful analysis of the coefficients Cn{k) in (^) shows, all but the first few are combinations 
of long orbits with short orbits approximating the long ones (Cvitanovic and Eckhardt 1989, 




(8) 
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Fig. 2: Scaling of tlie coefficients in the cycle expanded Selberg zeta product (H). 
Shown is log^g \Cn{k = 0)| vs. n for the same distance to radius ratios as in Figure 1. 



Artuso, Aurell and Cvitanovic 1989a,b, Eckhardt 1992, Eckhardt and Russberg 1992). If the 
index a^b describes the long orbit composed of n repetitions of a shorter one a and a second orbit 
b, then the expansion (P) will contain contributions of the form (e.g. for j = 0) 



a"b\ 





Aa'ife 


l/2\ 









(9) 



The second term in the parenthesis contains only terms comparing the long orbit with the two 
approximants. In hyperbolic systems, orbits are by definition exponentially unstable. Thus to 
stay near an orbit for n traversals requires one to be exponentially close to it. Thus the difference 
between orbits a"6 and a"~^6 is one more traversal very (exponentially) close to a. Thus the 
differences in action are exponentially close to zero, the ratios in instabilities close to 1. If now 
the Maslov phases (a discrete quantity) fit as well (one can assure this by a proper choice of 
symbolic organization of orbits, Eckhardt and Wintgen 1991, Eckhardt 1992), all that enters in 
(^) are the deviations from one, i.e. an exponentially small quantity. This is enough to turn the 
expanded product (|^) convergent for real energies. 

This is a ffist example of order among all this chaos: long periodic orbits are not random 
numbers but are strictly correlated with shorter orbits, so that one can actually estimate the 
contributions from long orbits rather accurately and obtain convergent results for the density of 
states and related quantities. 

For large n the decay is only exponential rather than superexponential, in contrast to the 
classical Selberg zeta function (Cvitanovic and Eckhardt 1991). Techniques have been developed 
to improve this further to recover faster than exponential convergence (Eckhardt and Russberg 
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1992). A new semiclassical Selberg type product, similar to the classical one, is currently under 
investigation. Preliminary test indicate improved convergence. The calculations shown here were 
performed with the expansion (Q) directly, using orbits of period up to 13. 



One test of the semiclassical expression is to compare zeros of (§) with poles of the S-matrix. 



The result (Cvitanovic and Eckhardt 1989, Eckhardt et al 1992) is that the differences between 
exact and semiclassical resonances decrease with increasing wavenumber, i.e. when approaching 
the semiclassical limit. This is in spirit with the semiclassical approximation and confirms that 
Gutzwillers (1990) analysis is a bona fide semiclassical theory. 

A plot of the time delay for a larger interval in k shows large fluctuations, and lots of irreg- 
ular features (Fig. 3). They come about because every resonance of the S'-matrix contributes a 
Lorentzian to the time delay. If the positions of the resonances are ki = Si — iji, then 



Cross sections differ from the above expression in that the weights of the Lorentzians are not 
unity (Ericson 1960, Brink and Stephen 1963) and also show fluctuations. Signals of this kind 
were extensively studied in nuclear physics (Ericson 1960, Brink and Stephen 1963) in connection 
with the compound nucleus. It was shown that a correlation function contains information on the 
average widths of resonances and thus on the lifetime of the compound nucleus. A semiclassical 
interpretation has been attempted by Bliimel and Smilansky (1988), starting directly from the 
S'-matrix. In this limit, however, the same information should also be contained in the time 
delay and its correlation function (for a discussion of deviations in extreme quantum cases, see 
Lewenkopf and Weidenmiiller 1991). 



Consider the fluctuating part of the time delay (after subtracting the mean), more precisely a 
segment Ts{k) = Tf{k)w{k), projected out with a window function w{k) (e.g. a box of width Ak 
centered around k^ or a Lorentzian). The autocorrelation function of this segment is then given 
by the integral 



{k-s,r + jf 



(10) 



CORRELATION FUNCTIONS 




(11) 



Its Fourier transform K{X) is the structure function. 




(12) 
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Fig. 3: Time delay for scattering off tfie three disk system at d : R = 2.5. Tlie lower 
frame shows a magnification and comparison between different maximal periods of 
orbits included. 
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Fig. 4: Autocorrelation function of tlie time delay for d : R = 2.5. It was computed for 
intervals of An = 50 starting at the fc-values indicated in the figure. For comparison, 
the bold line shows a Lorentzian of width 7 = 0.15 (Eq. p!^). 



Since the time delay is a superposition of Lorentzians, the correlation function for small k 
probes their widths. In a simplified analysis, assume that the window w{k) selects a finite number 
of resonances in each T{k). Then 

c{k) = yy^ fdk - — — . - — — . (13) 

ijJ {k- k/2- Sif + ^i^ {k + k/2 - s,y + ^ ' 

which, in the diagonal approximation, becomes 

^('^) = ^ (k - s^f + ^il + ■ ^^^^ 

The last form is an approximation by a single Lorentzian of effective width 27. 

Comparison with the computed autocorrelation function (Figure 4) shows that a Lorentzian 
describes the correlation function only for very small k. For larger ones, C{k) even becomes 
negative. A similar discrepancy was noted by Wardlaw and Jaworski (1989) in their analysis 
of the phase shift for a particle scattered off a leaky surface of constant negative curvature (a 
model introduced by Gutzwiller 1983). In a later study, Shushin and Wardlaw (1992) could 
actually improve on the form of C{k) by adding correlations between the positions of resonances, 
predicting a correlation function of the form 

gw(^)^7 %^,'~';' • (15) 



Fig. 5: Structure function for the time delay sliown in Figure 3. The inset shows the 
structure function over a larger range on a logarithmic scale. Clearly visible is the 
exponential fall off related to the Lorentzian of the autocorrelation function. 



This form does not fit the data any better, presumably because the underlying ensemble of 
resonances is different from the assumed GUE ensemble (see e.g. Haake 1991). 

The large k behaviour of the correlation function is difficult to access from the resonance 
representation (p!OD. This behaviour on the other hand is reflected in the small A behaviour of the 
structure function. Using the semiclassical representation (0) for the time delay, the structure 
function for a segment Ts{k) becomes 



K{\) 



V V ^(A - tL,) 

rhj\det{l-M;)\ 



EE 



^2 



r,2 



p r- 



\ I det(l - Ml 



(A - rLp) + ApAp'W{X - Lp)w{X - L 



pj ' 



(16) 
(17) 



p p' 



where w is the Fourier transform of the window function; P and P' denote periodic orbits (perhaps 
multiple traversals of shorter ones) with weights Ap, Api and of length Lp and Lpi, respectively. 
If the window is sufficiently wide, then the Fourier transform will be narrow and no terms will 
overlap when taking the absolute value squared. Therefore, only the diagonal contributions in 
([ITD survive and the structure function will show individual periodic orbits for small A. For 
very large A, interferences between different periodic orbits will eventually turn the sum into the 
exponential decay required by the Lorentzian result for the small r behaviour of the correlation 
function. 



RELATION TO BOUNDED SYSTEM 
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As the disks are moved close enough to touch, they enclose a region of the shape of a tipped 
triangle. Since this is a bounded system, it has a discrete spectrum. The Wigner time delay then 
turns into the density of states of the bounded system and the resonances become infinitely sharp. 
Semiclassical quantization attempts for this system have been reported elsewhere (Tanner et al 
1991). Here I would like to describe the relation between the correlation function and structure 
function for the time delay as found above to the semiclassical theory of two point correlation 
functions by Berry (1985). 

The small A behaviour of the structure function with the periodic orbits is in agreement with 
Berry's semiclassical analysis and numerical results (e.g. Wintgen 1987). The large A exponential 
decay is consistent with it, since the decay rate is related to the width of the resonances; in a 
bounded system, the widths goes to zero and K{X) settles to a constant (Levandier et al 1986, 
Pique et al 1987). What seems to be missing is the intermediate regime, where in a bounded 
system K{X) should increase linearly. 

For A large enough so that the periodic orbits overlap, but not too large so that the cancellations 
take over, the correlation function near A sums over all orbits with length near A, 

E I ,.,,5 ...m "^^"''- (18) 



near A 



det(l - M;; 



The approximate form is valid for intermediate A and contains the classical escape rate F (Kadanoff 
and Tang, 1984, Cvitanovic and Eckhardt 1991), up to a factor L. In the case of bounded sys- 
tems, F = and K{X) increases linearly with A. For open systems, it attains a maximum near 
A = 1/F. For the three disk system at d : R = 2.5 : 1, F ^ 0.7, so that it is impossible to 
detect the intermediate classical region. However, the exponent F clearly does not describe the 
exponential decay for large A, as speculated by Bliimel and Smilansky (1988). Incidently, notice 
that in this case the quantum decay is slower than the classical decay. Thus the interferences 



between off-diagonal terms in (17) act to increase C(A) compared to the classical value and not 



to decrease it as in the case of the bounded system (Berry 1985). 

CONCLUSIONS 

Investigations in 'quantum chaos' and in particular on Gutzwiller's periodic orbit formula for 
the density of states have helped to uncover a considerable amount of order among irregular 
looking quantum data. The main tool is the Fourier transform which facilitates the transition 
from the energy domain to the time domain, where the dynamics is dominated by correlations; 
they give rise to longer range correlations in energy. I have demonstrated these ideas using as an 
example a simple scattering system, where one has good control over both classical and quantum 
dynamics and where it is possible to sum the semiclassical expressions, despite the exponential 
proliferation of orbits. The results confirm that the cycle expansion of Selberg products (|^) 
provides a convenient and accurate method for the calculation of periodic orbit sums in hyperbolic 
systems. It has allowed for a accurate test of Gutzwiller' semiclassical theory. 
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A number of problems still remain. Among them are the extension of periodic orbit summa- 
tion techniques to non-hyperboUc cases, the summation of expressions for cross sections and the 
extension to more degrees of freedom than just two (perhaps even field theories). 1 am optimistic 
that eventually, a semiclassical theory for chaotic systems, perhaps as useful and reliable as the 
well-known WKB methods, will emerge. 
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